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Abstract. - Working together in groups may be beneficial if compared to isolated efforts. Yet 
this is true only if all group members contribute to the success. If not, group efforts may act 
detrimentally on the fitness of their members. Here we study the evolution of cooperation in public 
goods games on scale-free networks that are subject to deletion of links that are connected to the 
highest-degree individuals, i.e., on networks that are under attack. We focus on the case where 
all groups a player belongs to are considered for the determination of payoffs; the so-called multi- 
group public goods games. We find that the effect of link deletions on the evolution of cooperation 
is predominantly detrimental, although there exist regions of the multiplication factor where the 
existence of an "optimal" number of removed links for the deterioration of cooperation can also 
be demonstrated. The findings are explained by means of wealth distributions and analytical 
approximations, confirming that socially diverse states are crucial for the successful evolution of 
cooperation. 
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Introduction. — Understanding the evolution of co- 
operation is a fundamental problem in social and biolog- 
ical sciences [Tl[2]. As a common theoretical framework, 
the prisoner's dilemma game [3] has received ample at- 
tention in the study of cooperation between selfish indi- 
viduals via pairwise interactions. However, cooperation 
is observed not only in pairwise interactions, but indeed 
even more so in groups involving more than two individ- 
uals. Examples are many, including public transportation 
[4] and attempts of averting global warming [5]. Public 
goods games (PGGs) are established, both theoretically 
as well as experimentally, to capture the essence of the 
dilemma that underlies such cooperative behavior. In a 
typical PGG experiment, n players are asked to invest 
into a common pool. They all know that the accumulated 
contributions will be multiplied by a factor r > 1, and 
then equally divided among all participants irrespective 
of their investment. If all players cooperate they increase 



their initial input by (r — l)c, where c is the invested con- 
tribution, i.e., the cost of cooperation. However, every 
player is faced with the temptation to defect by exploit- 
ing the contributions of other players and withholding its 
own. In particular, if a player chooses to defect while oth- 
ers cooperate it will get a higher net payoff. It is thus 
more likely that this strategy will spread in an environ- 
ment governed by natural selection due to the higher fit- 
ness of defectors. The conflict escalates as more and more 
players adopt the defecting strategy. In doing so, they ne- 
glect collective interests, which ultimately leads to social 
downfall and poverty, known also as the "tragedy of the 
commons" |6]. Human experiments reveal, however, that 
the level of cooperation is surprisingly high in one round 
games, but deceases fast from round to round if a game is 
played repeatedly [7]. 

Thus, in PGGs individuals face the temptation to reap 
benefits on the expense of others, and in doing so likely in- 
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duce further defecting actions which harm collective well- 
being. In contrast, cooperators are able to resist these 
temptations with the aim of facilitating group benefits and 
social welfare. Defectors therefore bear no costs when col- 
lecting identical benefits as cooperators, which ultimately 
leads to widespread defection. This theoretical prediction, 
however, disagrees with experimental findings [S]. To re- 
solve this disagreement between observations and theory, 
several mechanisms have been proposed that promote co- 
operation. Hauert et al. [9] introduced voluntary partic- 
ipation in PGGs and found that it supports cooperation 
in a Red-Queen-like manner. Szabo et al. [TU] also stud- 
ied the impact of voluntary participation in PGGs on a 
square lattice, reporting that the introduction of loners 
leads to a cyclic dominance of the three strategies. Pun- 
ishment too has been identified as a viable route to co- 
operative behavior [TTHT^ . although its effectiveness de- 
pends on whether the participation in the PGG is optional 
[14] , and on whether the interactions amongst players are 
structured [15] . The reward as a means to cooperation in 
spatial public goods games has also been considered [16]. 
Recently, social diversity introduced via heterogeneous in- 
teraction networks and participation of players in multiple 
groups [T7l[T8] was identified as a viable route to cooper- 
ative behavior in PGGs, as was the relaxation of strategy 
adoption criteria to account for mutation and random ex- 
ploration of available strategies [19]. Effects of different 
temporal and spatial scales in evolutionary games have 
also been revealed as being potentially important [20H22] . 
Fowler and Christakis observed experimentally that coop- 
erative behavior can spread among people with personal 
ties [23]. However, quantitative research on the effects of 
altering the social ties in multi-group PGGs is more diffi- 
cult to come by [23]. 

In this paper, we therefore study the evolution of co- 
operation in PGGs on scale-free networks [25H27] where 
players can contribute and receive payoffs from all the 
groups with which they are affiliated. More precisely, in 
a multi-group PGG, every individual participates also in 
the groups centered around its direct neighbors, as was 
proposed recently in [17]. The topology of the interac- 
tion network also plays an important role in the evolution 
of cooperation (for a comprehensive review see j28j). In 
particular, it is well-known that heterogeneous interaction 
networks promote cooperation if the interactions are pair- 
wise [29l[30]. Since, in a scale- free network, there exists 
a small set of nodes with many more links attached to 
them than to other nodes forming the network, it is im- 
portant to assess the role of this highly connected nodes by 
the evolution of cooperation. We are thus led to consider 
multi-group PGGs on scale-free networks that are subject 
to deletion of links connected to the highest degree indi- 
viduals. 

We shall examine quantitatively the effects of network 
heterogeneity on cooperation in multi-group PGGs. In 
particular, we generate scale-free networks according to 
the method in |31j and then remove a number of links con- 
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Fig. 1; Fraction of the largest connected component / in de- 
pendence on the number of deleted links niL- The original 
network size here was A'' = 5000, and the results were obtained 
as averages over 1000 different initial network configurations. 

nected to the highest degree nodes. The deletion of such 
links is traditionally referred to as network attack 02] , and 
it was shown that the resilience to such actions is crucial 
for situations such as fast information transfer within the 
world- wide- web ^33j , uninterrupted supply with electricity 
[33], fast spread of epidemics and viral infections [55H55] . 
and cascade failure in networks [3p]. Several studies elab- 
orated analytically on the resilience of complex networks 
via the usage of percolation theory [40H43] . Interestingly, 
we find that the deletions mainly impair the evolution of 
cooperation in multi-group PGGs. However, the effect is 
monotonously negative only within a rather narrow re- 
gion of the multiplication factor, whereas outside this re- 
gion there exists an "optimal" number of deleted links at 
which the fraction of cooperators on the network is mini- 
mal. We examine the resulting distributions of wealth in 
terms of the dependence on the number of removed links 
in order to shed light on the reported results. We confirm 
that socially diverse states are crucial for the successful 
evolution of cooperation, as was argued recently in [17j . 

The remainder of this paper is organized as follows. In 
the next section we describe the employed multi-group 
PGG, and in Section HI we present numerical results. 
Section IV features an analytical treatment that provides 
further insights and explains our findings. Lastly, we sum- 
marize the main results and outline their potential impli- 
cations. 

Multi-group public goods games on scale-free 
networks. — In reality, interactions among individuals 
can often be described by complex networks. Here we 
focus on scale-free networks whereby each vertex corre- 
sponds to a particular player x. The generation of the 
network starts with two connected players, and subse- 
quently every new player is attached to two old players 
already present in the network, whereby the probability 
H that a new player will be connected to an old player 
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Fig. 2: Schematic presentation of a multi-group public goods 
game. Player x is not only a member of the group in which it 
is focal (the group is marked with a dark gray circle), but it is 
also a member of four other groups (one such group centered 
around its neighbor y is marked with a light gray circle.) 



X depends on its degree according to 11 — kx/^ky. 
This growth and preferential attachment scheme yields a 
network with an average degree n = (1/-/V) of four, 

and a power-law degree distribution with the degree expo- 
nent of 3 |3T]. Before the PGG starts, we delete a total of 
rriL links that are connected to the highest degree nodes. 
According to Albert et al. 32 , this corresponds to an at- 
tack, and we are interested in how different values of 
affect the evolution of cooperation in multi-group PGGs. 
As shown in Fig. [l] with increasing of rriL there exists 
a transition from an integrated network (there exists a 
path between any pair of nodes) to a fragmented network, 
where some groups of nodes become isolated and cannot 
be reach from the largest network component. The frac- 
tion / — S/N can then be defined, where S is the size of 
the largest connected component and N is the size of the 
original network (here we have N = 5000). Note that for 
rriL < 4000, S fa N so that the original network remains 
essentially intact. For > 6000, however, S N. At 
the transition point the value of is about 5000 if the 
size of the network is = 5000, as was the case in Fig. [T] 

Initially each player on site x is designated cither as 
a cooperator {sx — C) or as a defector [s^ = D) with 
equal probability. Cooperators contribute c (cost of coop- 
eration) to each public goods game, while defectors con- 
tribute nothing. Subsequently, the total contribution that 
accumulates within a group is multiplied by the multi- 
plication factor r > 1 that reflects synergistic effects of 
cooperation, and the resulting amount is divided equally 
among all the players irrespective of their strategies (and 
initial contributions) . 

Since players participate in multiple groups, each player 
X acquires its total payoff Px not only by participating in 
the group consisting of its direct neighbors, but also in 
all the groups that are centered around its direct neigh- 
bors (see Fig. TTl. Figured] demonstrates such a setup 
schematically, where player x is a member of five groups, 
having five members each. One is the group centered 
around player x (marked with a dark gray circle), but 
there are also four other groups centered around its neigh- 
boring players in which player x is also member. An ex- 
ample of such a group is marked with a light gray circle 
and is centered around player y. The payoff of player x 



with strategy Sx that is associated with the group centered 
around player y is then 
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(1) 



where z — stands for y, Sz is the strategy of the neighbor 
z of y, and k^ is its degree. The strategy is s^, = 1 if x 
cooperates and Sx = if it defects. The cost c refers to 
the contribution of every cooperator in each group it is a 
member. For example, if a cooperator is member in five 
groups, then its total cost will be 5c, i.e., one c for each 
group (or equivalently, one c for each public goods game 
played. Without loss of generality, the cost c is set equal to 
1. The total payoff of player x comes from the summation 
over all the groups 
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where r^, denotes the groups that are centered around the 
neighbors y and the player x itself. The number of groups 
in which player x participates is thus kx + 1. 

At each time step each player x acquires its payoff 
Px as described above. Subsequently, all players update 
their strategy synchronously by selecting at random one 
of their direct neighbors y. li Px > Py player x keeps its 
strategy Sx, but if Px < Py player x adopts the strategy 
of player y with the probability 



W(Sy Sx) 



P — P 

^ y ^ X 

AP 



(3) 



where AP^ax ensures the proper normalization and is 
given by the maximal possible payoff difference between 
players x and y [T7] . Results presented below are obtained 
on networks hosting N = 5000 players. Equilibrium frac- 
tions of cooperators pc are determined within 5000 time 
steps after sufficiently long transients were discarded. 

Simulation results. — First, we investigate the de- 
pendence of the fraction of cooperators pc on the enhance- 
ment factor r for different numbers of removed links at- 
tached to the main hubs to^, as shown in Fig. [3ja). It 
can be observed that the evolution of cooperation due to 
link deletions is affected strongly in multi-group PGGs. 
In particular, the cooperation at a given value of r as well 
as the emergence of complete cooperator dominance both 
vary significantly in their dependence on rriL- Indeed, if 
rriL — 4000 links are removed from the network before the 
multi-group PGG starts, the evolution of cooperation will 
be greatly impaired practically across the whole span of r, 
especially if compared to the niL — case. We can thus 
conclude that the decrease in network heterogeneity due 
to link deletions acts predominantly detrimental on the 
evolution of cooperation in multi-group PGGs. 

To examine the impact of link-targeted attack on the 
evolution of cooperation in PGGs more precisely, we show 
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Fig. 3: (a) Fraction of cooperators pc in dependence on the 
multiplication factor r for different numbers rriL oi deleted links 
connected to the highest degree players, (b) Fraction of coop- 
erators pc in dependence on m_t for different multiplication 
factors r. Note that for certain values of r there exists an "op- 
timal" number of deleted links rriL = fn'^^, at which defection 
thrives best. This effect, however, is quite subtle (see main 
text for details). The vertical Eixis has a logarithmic scale. 



in Fig. [3l^b) the dependence of pc on rTi^ for different val- 
ues of r. It can be observed that Hnk deletions have a 
predominantly negative impact on the evolution of coop- 
eration. More specifically, for low values of r the decrease 
of cooperation is most severe, which is in agreement with 
the fact that low multiplication factors are most challeng- 
ing for the evolution of cooperation, and accordingly, the 
decrease in the network heterogeneity due to link deletions 
is there most notable. It can also be observed that for low 
values of r there in fact exists an "optimal" number of 
deleted links at which defection thrives best (cooperators 
are most decimated). On the other hand, for larger r the 
effect of link deletions is monotonously negative, i.e., the 
more links are deleted the smaller the level of cooperation 
on the network. For r exceeding a threshold (r > 3.74 
|44j , which is linked to the survivability of cooperators on 
the square lattice), however, the negative impact of link 
deletions practically vanishes. This is because at such high 
values of r the synergetic effects of cooperation suffice to 
withstand defector attacks even in the complete absence 
of network heterogeneity. 

The reported decrease of cooperation upon link-targeted 
attack is consistent with previous results in that the con- 
nections among hubs facilitate cooperation in two-player 
networked games, such as the prisoner's dilemma or the 
snowdrift game [29l|45], as well as in multi-player PGGs 
[17j . In particular, a defector hub is frequently weak due 
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Fig. 4: The existence and dependence of the "optimal" number 
of deleted links m^^*, where defection thrives best, in depen- 
dence on the multiplication factor r. There exists only a rather 
narrow region (hatched gray) where the impact of link dele- 
tions is monotonously negative (m^^' does not exist). Outside 
this region, however, there always exists m°j^^ at which pc is 
minimal. 



to a negative feedback effect caused by its defecting neigh- 
borhood. Since there is then nobody left to exploit, a co- 
operator linked to such a defecting hub can overtake it 
and thus recover cooperation in a notable portion of the 
network. 

Turning to results presented in Fig. |3l[b) , we elaborate 
further on the existence of an optimal number of link dele- 
tions, denoted here by m"]^* , at which defectors are most 
widespread. Results presented in Fig.|3]indicate that there 
exists only a rather narrow range of r values (hatched 
gray), at which the link deletions have a monotonously 
negative impact on the evolution of cooperation. Out- 
side this region, however, there always exists a non-zero 
value of m'^^ at which pc is minimal. We note, neverthe- 
less, that following what is usually a sharp decrease of pc 
[especially for low values of r; see Fig. [3ljb)], the subse- 
quent increase of pc upon increasing rriL past m'^* is very 
modest indeed. The effect is mainly due to the creation 
of very isolated (although not fully disconnected) groups 
of cooperators that can prevail against defectors due to 
their isolation. This can be interpreted as group selection, 
although in our case the creation of groups is a sole con- 
sequence of the sparse network structure that emergences 
when TO^ approaches the transition point for network dis- 
integration (see also results presented in Fig.[T]). 

Finally, to gain further insights, we examine the result- 
ing wealth distributions for different values of tol and r, 
as presented in Fig. El It can be observed that multi-group 
PGGs induce heterogeneous wealth distributions that are 
characterized by a fat tail. The algebraic wealth distribu- 
tion agrees well with previous results pTlfTSlHS] . Notably, 
as the number of removed links increases, there emerges 
a cutoff in the distribution for the largest total wealth. 
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Fig. 5: Wealth distributions resulting from multi-group PGGs 
for different numbers of deleted links rriL connected to the high- 
est degree players. Values of the multiplication factor are: (a) 
r = 2.5, (b) r = 3.0 and (c) r = 4.0. 



The resulting deviation from the scale- free distribution ob- 
tained for rriL — induces a downfall of cooperation, espe- 
cially for small values of r, as evidenced in Fig. [2Ub). It is 
next of interest to summon our findings reported thus far 
and provide an analytical treatment to corroborate them. 

Analytical treatment. — In what follows, we pre- 
dict the payoffs of agents with respect to their degrees as 
well as the corresponding wealth distribution by virtue of 
an analytical approach. Assuming that cooperators are 
randomly distributed on the network and the fraction of 
cooperators is given by pc, we can express the payoff pij- 
as 



PC 

h + l 



(4) 



where Aji is the adjacency matrix of the underlying net- 
work. Neglecting the degree-degree correlation, we have 



N 



^ h + l ^ ^ k' + 1 ^ ^ 
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k'P{k') 1 
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(5) 



where (•) denotes the average of all nodes, P{k') is the 
degree distribution of network nodes, P{k'\kj) is the joint 
degree distribution defined by P{k'\kj) — k'P{k')/{k), and 
the identity YltZZ P{k')k' /(k'+l) = {k/{k + 1)) has been 
used. Substituting Eq. ([5]) in Eq. (|4|), we obtain 
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The total payoff is then given by 
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For fci 3> 1, we can obtain 
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Here we can have A> B, B > —1 and fc^ ^ 1, is thus 
proportional to ki. Distribution of wealth P{P) can then 
be obtained by P{k)dk = P{P)dP for P{k) = {kf/2k-^: 



P{P) 



^^/, s dk 
^ ^dP 
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(10) 



Since in the mean-field approximation the payoffs of 
hubs are proportional to their degrees, the connections 
among hubs play a key role in the maintenance of coop- 
eration pn[29ll45] . The deletion of links among hubs thus 
prevents cooperation, as shown in Fig.[3][b). 

Summary. — We have investigated the evolution of 
cooperation in multi-group PGGs on scale-free networks 
that were subject to link-targeted intentional attack. The 
main finding is that the deletion of links connected to 
hubs has a predominantly negative impact on the evo- 
lution of cooperation when players can participate in mul- 
tiple groups. However, the effect is strictly negative (the 
more links are removed the lower the fraction of coopera- 
tors) only in a rather narrow region of the multiplication 
factor, whereas outside this region there exists an interme- 
diate number of removed links at which defection thrives 
best. The latter effect, however, is quite subtle and is due 
to the emergence of group selection that sets in because of 
the sparse network structure that occurs when large num- 
bers of links are removed. The results can be explained by 
means of wealth distributions and an analytical treatment 
that is further supported by heuristic arguments, albeit 
the subtleties of the existence of an intermediate number 
of link removals for highest levels of defection cannot be 
precisely accounted for (due to the emergence of a rather 
contrived version of group selection, as described) . To sum 
up, we confirm that socially diverse states are crucial for 
the successful evolution of cooperation, where the upper 
limit of personal wealth introduces a cutoff that warrants 
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the most challenging environment for the spread of de- 
fection. Given these facts, strong parallels can be drawn 
between the evolution of cooperation on complex networks 
and processes such as spread of viral diseases or the effec- 
tiveness of information retrieval, which were also found to 
be strongly affected by intentional attack [35]. We hope 
that our work will inspire further research in this direc- 
tion and that it will serve as a useful source of information 
when striving towards the understanding of the stability 
of cooperation under attack in social networks. 
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